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Relativistic dissipative hydrodynamic equations are extended by taking into ac- 
count particle number changing processes in a gluon system, which expands in one 
dimension boost-invariantly. Chemical equilibration is treated by a rate equation 
for the particle number density based on Boltzmann equation and Grad's ansatz for 
the off-equilibrium particle phase space distribution. We find that not only the par- 
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tide production, but also the temperature and the momentum spectra of the gluon 
system, obtained from the hydrodynamic calculations, are sensitive to the rates of 
particle number changing processes. Comparisons of the hydrodynamic calculations 
with the transport ones employing the parton cascade BAMPS show the inaccuracy 
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of the rate equation at large shear viscosity to entropy density ratio. To improve the 
rate equation, the Grad's ansatz has to be modified beyond the second moments in 
momentum. 

PACS numbers: 47.75. +f, 24.10.Lx, 24.10.Nz, 12.38.Mh, 25.75.-q 



* el@th. physik. uni-frankfurt.dc 

t Azwinndini. Muronga@uct.ac.za 

* xu@th.physik.uni-frankfurt.de 

§ Carsten.Greiner@th.physik.uni-frankfurt.de 



I. INTRODUCTION 



Relativistic dissipative hydrodynamics has established itself as an effective theory for 
investigations of phenomena in relativistic heavy-ion collisions owing to its success in de- 
scribing the elliptic flow Vi data measured at the Relativistic Heavy Ion Collider (RHIC) at 
Brookhaven National Laboratory (BNL) [lHa] • Most hydrodynamic approaches [2tJ, IT], [8( do 
assume an instantaneous chemical equilibrium of matter constituents. However, this requires 
an infinite large transition rate of microscopic processes which drive the system towards the 
chemical equilibrium and thus breaks the consistency of the physical description for matter 
with a finite viscosity. Another extreme but physically consistent case, as considered in 
Ref. 9J, is to keep the number of the constituents constant by assuming that there are no 
particle number changing processes. In this case deviation of the one-particle phase space 
distribution from its chemical equilibrium form becomes larger and larger when the system 
expands. In real ultrarelativistic heavy-ion collisions bremsstrahlung processes of gluons 
and quarks, also gluon fusion to quark and antiquark, and the reverse processes govern the 
chemical equilibration of gluons and quarks 1CH14J|. Therefore, particle number changing 
processes with realistic rates should be included in underlying hydrodynamic descriptions 
for quark-gluon plasma (QGP). 

Because the energy density e ~ AT 4 , where A is the fugacity and T the temperature, for 
keeping the energy density unchanged the temperature decreases by 20% when the fugacity 
increases from 0.5 to 1. Thus, an effect of chemical equilibration is expected to be pre sent 
for the temperature, which is an important detail, for instance for the dilepton ll|, ll5Ml7| 
and photon 18j yields in ultrarelativistic heavy-ion collisions. Moreover, to what degree 
chemical equilibration of quarks and gluons is achieved at the phase transition might also 



be crucial for modelling hadronization using recombination models 19]. 

In this work we develop a dissipative hydrodynamic approach including particle number 
changing processes via a rate equation to investigate the bulk properties of the QGP at 
RHIC As a first step we simplify the expansion dynamics at RHIC by assuming a one- 
dimensional expansion in the beam axis with Bjorken boost invariance. The QGP is con- 
sidered to possess a constant rj/s value during its evolution. Different from the previous 



studies of chemical equilibration of the QGP basing on either ideal 
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hydrodynamic equations, we apply the second-order Israel- Stewart 
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11 ] or first-order 



as well as the 



extended third-order hydrodynamic equations [22] to describe the QGP evolution. The hy- 
drodynamic equations are now coupled with a rate equation for the particle number density. 
The microscopic interactions, which are responsible for the viscous behaviour and change 
of the particle number, are simplified by elastic binary 2 — > 2 and inelastic multiplication 
and annihilation 2 -H- 3 processes with isotropically distributed collision angles. Although 
the evolution of multi-component systems as discussed in Ref. 23] is gaining an increasing 
attention, here we consider only a pure gluon system for the sake of simplicity. 

To examine the applicability of the new description we compare the hydrodynamic solu- 



tions with those calculated from a parton cascade, the Boltzmann approach o 



scatterings (BAMPS) [13 
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in a similar manner as done in Refs. [9|, 
tj/s is extracted from the BAMPS calculations using the procedure introduced by us in Ref. 

21]. 

The paper is organised as follows. In Sec. [IT] we introduce the dissipative hydrodynamic 

equations for a one-dimensional boost-invariant expanding system of gluons in presence of 
particle production and annihilation processes. Solutions are compared with those without 
the inelastic processes and those assuming chemical equilibrium. This quantitatively demon- 
strates the effects of chemical equilibration on physical observables. Influences of the initial 
conditions on these observables is studied in Sec. IIHI In Sec. |IV]we compare the results of 
hydrodynamic calculations with those from BAMPS. Conclusions are given in Sec. |VJ 

II. VISCOUS HYDRODYNAMIC EQUATIONS AND RATE EQUATION 
DESCRIBING CHEMICAL EQUILIBRATION 



Derivations of dissipative hydrodynamic equations either from the entropy production 



principle or the moments method were reported in Refs. 22|, |26|, |27(. Ignoring the bulk 



viscosity and the heat conductivity, one obtains the evolution equation for the shear tensor 
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and A Q/ g = g a/3 — u a up with the metric g a p = diag(l, —1, —1, —1). u^ is the fluid velocity 
with u^u^ = 1. r n = 2r//3 2 is the relaxation time, where r\ denotes the shear viscosity and 



02 — 9/(4e). The "IS" part is exactly the equation from the Israel- Stewart theory, while 



the "third-order" part indicates a higher order correction term derived in Ref. 22j. We 
obtained a = —8/9. 

We simplify the dynamical evolution of a QGP created in ultrarelativistic heavy-ion 
collisions by a one-dimensional boost-invariant expansion of massless Boltzmann particles. 



For this case w M = (t, 0, 0, z)/t where r = \Jt 2 — z 2 is the proper time. With d^u^ = 1/r 
and u^d^ = 4- Eq. (pQ) reduces to 

n 8 e l7f it 1 e 3tt 2 3 vr 2 T 3 vr 2 e vf 2 

vr = 1 h -vf— + -7f- H 1 1 4 — , (3) 

TV 27 r 2r 2 T 2e 2 ex 2eT 2ee er ' v; 

where fr = tt 33 = — 2tt 1:l = — 2tt 22 is the only independent component of i\^ v in the local 
rest frame. The dot denotes the derivative with respect to r. Time evolution equation for 
the energy density e follows from the conservation of energy, d^T^ = 0, where T^ u is the 



energy-momentum tensor. For a one-dimensional boost-invariant system one obtains 221.127] 
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Equation (E} corresponds to the Grad's ansatz [22] for an off-equilibrium distribution of 
particles in phase space 

3tt (\ 



f(x,p) = f eq (x } p)X 
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where f eg (x,p) = ge~ E ' T is the equilibrium distribution, g is the degeneracy factor and g = 
16 for gluons. A denotes the fugacity and quantifies the degree of the chemical equilibration. 
The distribution (J3J) satisfies the matching of the particle number and energy densities, 
n = Xn eq = \gT 3 /ir 2 and e = Xe P q = 3\gT 4 /ir 2 , by which an off-equilibrium state is matched 



to a fictitious equilibrium state 28|. The matching conditions allow to define temperature 



for an off-equilibrium system. Assuming chemical equilibrium, i.e., A = 1, one obtains 

With this one can solve the coupled equations (E]) and (J3D once the initial conditions are given. 
The assumption of chemical equilibrium is made in most of the hydrodynamic calculations 
[2J-I4J, |7|, ISJ . However, an instantaneous and complete chemical equilibrium requires an infinite 
transition rate of particle number changing processes, which is not consistent with a finite 



shear viscosity. Without the assumption of the chemical equilibrium, the temperature 
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depends also on the particle number density. Thus, we need an additional equation for the 
time evolution of n. 

We assume that the space-time evolution of gluons obeys the Boltzmann equation 



P< 1 d fl f(x,p) = C[f] 
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where C[f] is the collision integral and contains microscopic details of particle interactions. 
Comparisons between the solutions of the hydrodynamic equations and those of the Boltz- 
mann equation will be shown in Sec. IIVI Here we use the Boltzmann equation to derive an 
equation for the particle number density n. 

Integration of Eq. (|SJ) over dT = d 3 p/[(2ir) 3 E] leads to 



dTp^f(x,p) = f dTC[f]. 



(9) 



The left hand side can be rewritten as the derivative of the particle number current, d^N 11 , 
which is equal to h + n/r for a one-dimensional boost-invariant expansion. The right hand 
side presents the source of particle production and annihilation. Considering two to three 
particles and vice versa as the only particle number changing processes, one obtains the rate 
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where -R23 and R32 denote the collision rates per particle for inelastic 2—7-3 and 3—7-2 



processes. The 



are given by 30 



actor 1/2 and 1/3 indicate that colliding particles are identical. The rates 
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Vrei — (j?i +P2) 2 : /(2EiE 2 ) is the relative velocity of two colliding particles with four momenta 
P\ and p 2 . The general definitions of the cross section a 23 and I 32 can be found in Refs. 
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3l| . For interactions with isotropic distributions of collision angles I 32 is related to a 23 



via I 32 = 1927r 2 cr 2 3 |l3j. Using the approximate distribution (J5J) we obtain 

1 „ 1 „ 1 



which leads to 



^nR 23 - ~nR 32 = ^ 2 (1 - \)a 23 , (14) 
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Equation (ITS"]) is similar to the rate equations derived in Ref. 10] . 



The hydrodynamic equations (J3J), (05, and (ITS"]) describe kinetic and chemical equilibration 
of a gluon matter. Whereas the value of rj does not depend on the types of microscopic 
interactions [30j, the net particle production, which explicitly affects the shear pressure ft 
and n, is of course strongly dependent on the strength of inelastic collisions. 

The standard viscous hydrodynamic approaches like the ones employed in Refs. [3|, U, |2j 
are up to second order in gradients and dissipative quantities. In the one-dimensional case 
studied here, the relevant gradient is the expansion scalar <9 M u M = 1/r and the relevant 
dissipative quantity is the shear pressure ft. Both r n /r and ft/e must be small to ensure the 
validity of the hydrodynamic approach. In our notation a term of g-th order has the form 
(r 7r /r)' 3l (7f/e)' 32 with q 1 + q 2 = q. To make a proper order counting for Eq. ([3j) we rewrite 
T/T = e/e — n/n according to Eq. (jj|) and insert Eqs. fl3| and (fl5|) for e and h into (]3]). 
We then multiply Eq. (]3]) by r^/e. Keeping terms up to second order we obtain 

7f 4 7f 8 e 1 _ 
- = ---- + ---^(1-^3. (16) 

We refer to Eq. ffTB"]) as the Israel-Stewart (IS) equation. A similar equation, used in Refs. 
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9l. |29| assuming particle number conservation, i.e., a 23 = 0, contains the term 7f 2 /(er) which 
is neglected here, because this term times T n /e is of third order. Keeping terms up to third 
order we obtain 

ft Aft 8 e „7f 2 1 / „7f\ ,_ ,. . „. 

We refer to Eq. ( IT7|) as the "third-order" equation. 

In the following we solve the IS equation ( TIB"]) as well as the third-order equation ( 1T7|) 
coupled with Eqs. (]4]) and f lT5|) . We consider a gluon system with a constant 77/5 = 0.35, 
which is a rough estimate of the upper bound of the rj/s value found at RHIC [4j. The 



entropy density s is taken at kinetic equilibrium, i.e., at n = 

s = 4n-n\n\. (18) 

The fugacity is related to e and n via A = n/n eq = (tc 2 / g)nT~ 3 = (n 2 / g)n(e/3n)~ 3 . Since 
i] ~ e\ m f p , where \ m f P is the gluon mean free path, we have r]/s ~ (T/s)(n\ m f p ). With the 
choice 

a23 = ^7 = A(4-lnA)T* (19) 

and the same cross section of binary collisions o 22 = cr 23 , a constant rj/s = 0.35 can be 
obtained in kinetic transport calculations, as will be shown in FigJH in Sec. IIVI 

To demonstrate effects of the chemical equilibration on observables we compare results 
of calculations with non- vanishing a 23 to those with a 23 = (particle number conservation) 
and with A = 1 (instantaneous chemical equilibrium). In the case of a 23 = Eq. (jl~5|) gives 
h = —n/r and leads to n{r) = n(T )T /r. For A = 1 we have n = n eq = (3g/7r 2 ) 1//4 e 3 / 4 /3, 
which is the solution of Eq. (j!5p only if n = 0, i.e., rj = 0. This corresponds to infinite 
transition rates of particle number changing processes, which keep the system in chemical 
equilibrium. Here we also see that the setup of A = 1 and a non-vanishing shear viscosity rj 
is physically inconsistent. 

Figured] shows the results for a thermal initial condition at tq = 0.4 fm/c with T(tq) = 0.5 
GeV and A(r ) = 1. To see the effects of the finite shear viscosity and the particle number 
changing processes, we rescale the particle number density, the temperature, the energy and 
the entropy density by their time evolution in an ideal fluid. We first compare the results 
for non- vanishing a 23 (solid curves), a 23 = (dashed curves) and A = 1 (dotted curves) 
using the IS equations. With the particle number conservation (023 = 0) n ■ t/tq = n(ro) is 
constant. The non-vanishing a 23 leads to a net increase of particle number, because the finite 
shear viscosity brings the system away from chemical equilibrium and the system becomes 
undersaturated (A < 1). This is clearly demonstrated in the time evolution of the fugacity 
in Fig. [T] (b). It follows from Eq. f[T5"j) that the larger the value of a 23 , the larger is the 
particle number increase and thus the faster is the restoration of the chemical equilibrium. 
In the limit of instantaneous restoration (A = 1), corresponding to a 2 3 — > 00, a maximum of 
particle productions is achieved, as shown by the dotted curve in Fig. [1] (a). 

Due to the viscous effect, the decrease of the energy density is slower than in the ideal 
fluid. Three curves in Fig. [T] (c) differ only marginally. Same is observed for the ft/e ratio 
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IS, r|/s=0.35, o 23 =(X.(4-lnX.)T' i )" 
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FIG. 1: (Color online) Time evolution of (a) the rescaled particle number density n ■ (t/tq), (b) 
the fugacity A, (c) the rescaled energy density e • (t/to) 4 ' 3 , (d) the ratio 7f/e, (e) the rescaled 
temperature T ■ (t/tq) 1 ' 3 and (f) the rescaled entropy density s- (t/tq). The cases 023 = (dashed 
curves) and A = 1 (dotted curves) are calculated using the IS equations. The standard case with 
Eq. ([T2J) is calculated using both the IS (solid curves) and the third-order equations (dash-dotted 



curves . 



in Fig. [T] (d). Time evolution of the energy density and the shear pressure depend mainly 
on the viscosity and not much on the details of microscopic interactions. Formally, the weak 
dependence of e and n on (T23 or A is explained by the fact that the corresponding term in 
Eq. 02]) is proportional to T/T, which has a logarithmic dependence on the temperature. 

Because of T = e/(3n) the rates of the particle number changing processes influence 
the value of the temperature. The larger the rates, which are infinity in the case of the 
instantaneous equilibration (A = 1), the smaller is the temperature, as shown in Fig. He). 
Between the two limits, (J23 = and A = 1, there is a 50% difference at the final time r = 10 



fm/c. 

Due to the matching n = Xn eq and e = \e eq we obtain n ~ e 3 / 4 A 1 ^ 4 and thus s ~ 
e 3//4 A 1//4 (4 — In A) according to Eq. f lT8|) . Because the energy density is rather insensitive to 
the evolution of fugacity A [Fig. [I] (c)] and the function A 1/,4 (4 — In A) is almost constant 
for not very small A, the entropy density weakly depends on the fugacity, see Fig. Q] (f). A 
larger effect can be observed at late times in the case of (X23 = 0, because at that time the 
fugacity is quite small if particle number is conserved. 

The third-order correction terms to the IS equations reduce the entropy production. This 
is the main finding of Ref . 22] . In Fig. [1] (f ) we observe that the entropy density calculated 
from the third-order equations (dash-dotted curve) is considerably smaller than that from 
the IS equations (solid curve). Thus, in oder to obtain the same rj/s, a smaller shear viscosity 
77 is needed, which leads to smaller energy density and shear pressure, compared with the 
results of the IS equations. However, the interplay between these quantities and n, T, and 
A appears to be non-trivial. 
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FIG. 2: (Color online) Transverse momentum spectra according to the distribution ([5]). Depicted 
are the initial spectrum at tq = 0.4 fm/c and spectra at r = 1 fm/c, 4 fm/c and 8 fm/c. 

Since the rates of the particle number changing processes influence the temperature, they 
should as well influence the particle momentum spectrum, which is more closely related to 
experimental data. Using the distribution (jSJ) we calculate the particle transverse momentum 



10 

spectra at different times and demonstrate them in Fig. [2J Values of e, 7f, T, and A are 
taken from solutions of the IS equations in the three cases shown in Fig. [TJ According to 
Eq. 05]) the spectrum slope is approximately proportional to the temperature regardless of 
the dissipative correction. Significant differences in the spectra for the three cases at a late 
time, r = 4 fm/c or r = 8 fm/c, reflect the differences in temperature shown in Fig. [I] (e). 

III. INFLUENCE OF INITIAL CONDITIONS 

In ultrarelativistic heavy-ion collisions initially produced quarks are much less abundand 
than gluons. Chemical equilibration of quarks and gluons start with different initial values of 



fugacity. The initial fugacity of quarks is expected to be much smaller than one [lOj. In this 
section we repeat calculations performed for Fig. [I] but with an initial fugacity A(r ) = 0.2 
instead of 1. The initial condition is a kinetically equilibrated, but chemically disequili- 
brated system. The initial energy density is e(r ) = A(r )e e(? (ro) = (3g/7r 2 )A(To)[T(T )] 4 
with T(t ) = 0.5/[A(r )] 1/4 GeV. We calculate the IS equations, Eqs. (USD, ©, and (fT5l). 
for r]/s = 0.35. 

The results are shown in Fig. [3] for three cases: non- vanishing <7 2 3 [Eq. fTXTJT) ] by the solid 
curves, 023 = by the dashed curves, and A = 1 by the dotted curves. In the case of A = 1 
the initial temperature is 0.5 GeV; it is larger than 0.5 GeV in case Ao = 0.2 to hold the 
same initial energy density. The initial particle number density and the entropy density are 
larger in the A = 1 case compared with the values of the other two cases. 

As already observed in Fig. [H the energy density and the shear pressure, shown in Fig. 
[3] (c) and (d), are weakly affected by the value of 023. On the contrary, without particle 
number changing processes, (T23 = 0, the fugacity [see FigfS] (b)] decreases and the system 
goes far away from the chemical equilibrium, while the fugacity increases to 1 for large 023- 
Accordingly, pronounced differences are seen in the particle number density [Fig. |3](a)], the 
temperature [Fig. 12(e)], and the entropy density [Fig. |3](f)] between the non- vanishing cr 2 3 
and (T23 = cases. We also observe that the results of n, A and T in the non- vanishing (T23 
case relax to those of the A = 1 case. 

Comparisons in Fig. [3] underline the importance of microscopic details in applying hy- 
drodynamic approaches to describe both chemical and kinetic equilibration of an initially 
undersaturated system. 
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FIG. 3: (Color online) Same as Fig. Q] using the IS equation (fT6]) with kinetically equilibrated and 
chemically disequilibrated initial condition, A(ro) = 0.2. 

IV. COMPARISONS BETWEEN THE HYDRODYNAMIC AND TRANSPORT 

CALCULATIONS 



Kinetic transport theory describes the space-time evolution of particles by means of 
detailed treatment of microscopic processes. The dynamical behaviour of a particle system 
resembles the one described by viscous hydrodynamics if the particle mean free path is 
much smaller than macroscopic scales. Comparisons between hydrodynamic and transport 



calculations are a 
theories bl bl 22J 



ways useful to examine the applicability limits of viscous hydrodynamic 
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In this section we compare results from the new viscous hydrodynamic description in- 
cluding particle number changing processes presented in Sec. [TT] with results of transport 
calculations using the parton cascade BAMPS [l3[]. BAMPS simulates particle interactions, 
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especially inelastic 2 -h- 3 processes with full detailed balance, by using the stochastic inter- 
pretation of transition rates. 

From the kinetic theory point of view, the physical origin of a finite shear viscosity is the 
finite rates of microscopic collision processes. The relation between the (transport) collision 



rate and the shear viscosity coefficient was reported in Refs. 
formula derived in 
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, |32j. Here we use the 
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with P^ = JdTpi*p v C[f] and C = 3/(8eT 2 ). The collision integral C[f] contains details 
of microscopic interactions. We consider elastic binary 2 — y 2 and inelastic 2 -H- 3 processes 
with isotropic distributions of collision angles. The cross sections are parametrized by 

k 

° 22 = a23 = A(4-lnA)T 2 ' (21) 

where k is a parameter controlling the value of 77. With this choice one expects to obtain 
a constant rj/s as discussed around Eq. ( I19p . The transition probabilities of a 2 — y 2 
or 2 —y 3 process within a unit volume Av and a unit time step At are same and given 
by -P23 — -P22 — v re iCT23At J Av . The transition probability of a 3 — y 2 process is P 32 = 
(l/8E 1 E 2 E 3 )I 32 At/Av 2 with J 32 = 192vr 2 cr2 3 |13j. Figure H shows the 77/s ratio extracted 
from BAMPS calculations with k = 0.5, 1, 3, and 6 via Eq. ([20]) . rj/s decreases with 
increasing k. Each ratio is in good agreement with a constant. These constant rj/s values 
and corresponding a 23 parametrizations are used for solving the hydrodynamic equations 
(USD (IS) and (dU) (third-order) together with Eqs. © and flU*}. 

Figures El El [TJ El and |9] show the comparisons between the hydrodynamic and transport 
calculations for a one-dimensional boost-invariant expansion of a gluon system with the 
thermal initial condition T(tq) = 0.5 GeV, A(ro) = 1 at tq = 0.4 fm/c. The results are 
rescaled in the same manner as previously done in Fig. [H 

In Figs. El and El we observe that the energy density and the shear pressure from the 
third-order hydrodynamic calculations agree perfectly with the BAMPS results even for 
large rj/s values. In contrast, the viscous effect is overestimated in the IS equations. This 
becomes significant for large rj/s values, which mark the applicability boundary of the IS 
approach. These findings are in line with those in our previous work 22], where particle 
changing processes were not taken into account. 
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FIG. 4: (Color online) rj/s ratio extracted from BAMPS calculations with time-dependent cross 
sections. 
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FIG. 5: (Color online) Rescaled energy density from BAMPS, IS and third-order hydrodynamic 
calculations. 

The difference in the particle number density (Fig. [7]) between the hydrodynamic and 
BAMPS results is large at large 7]/s values. Although both IS and third-order equations 
give larger densities than those in BAMPS, the third-order results are closer to the ones 
from BAMPS . On the other hand, looking at the temperature T (Fig. [H]) and the fugacity 
A (Fig. [9]) we observe that the IS results show better agreement with the BAMPS results 
than the third-order ones. However, it is difficult to make conclusions about applicability of 
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FIG. 6: (Color online) Shear pressure to energy density ratio from BAMPS, IS and third-order 
hydrodynamic calculations. 
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FIG. 7: (Color online) Rescaled particle number density from BAMPS, IS and third-order hydro- 
dynamic calculations. 

a hydrodynamic approach basing on the observables T and A. These quantities are defined 
via e and n, but not solved directly from the hydrodynamic equations. If we consider the 
particle number conservation, the results on T and A from the third-order calculations are 
in very good agreement with those from BAMPS, as can be cocluded from our observations 
in Ref. 22]. In the situation considered here, small differences between hydrodynamic and 
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FIG. 8: (Color online) Rescaled effective temperature T = e/(3n) from BAMPS, IS and third-order 
hydrodynamic calculations. 
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FIG. 9: (Color online) Fugacity from BAMPS, IS and third-order hydrodynamic calculations. 

kinetic transport results on both key observables e and n translate into differences in T and 
A in a non-trivial way. 

To understand the differences in particle number densities between the hydrodynamic and 
BAMPS results, we examine Eq. (I14p . which is valid by virtue of the approximation fl5]). For 
this we calculate the left and the right hand sides of Eq. ( I14p using the actual values of R23, 
-R32, n, A, and 023 extracted from the BAMPS calculations. The results are shown in Fig. 
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FIG. 10: (Color online) Left and right hand sides of Eq. (|14p (denoted by LHS and RHS, respec- 
tively) calculated using the actual values of R23, R32, n, A and 023 from BAMPS. 

[TU1 Except for the case of 77/s = 0.05, n 2 (l — A)<723/2 is always larger than ni? 23 /2 — nR^/^, 
which leads to a stronger particle production in the hydrodynamic than in the transport 
approach, as seen in Fig. [JJ Therefore, it is obvious that the approximate distribution (j^J) 
must deviate from the one extracted from BAMPS. This deviation is demonstrated in Fig. 
[H] by the ratio of pr spectra at midrapidity extracted from BAMPS to the ones calculated 
using the Grad's approximation fl5J). The results are presented for rj/s = 0.35 at various 
times. The deviations of the ratios from 1 are on a level of 10% — 20% and become more 
pronounced at later times, r ~ 1 — 2 fm/c, at which the viscous effect is strongest, see Fig. 
EJ We thus need a modification of the Grad's approximation, which is a highly interesting 
project for future studies. 

V. SUMMARY AND OUTLOOK 



In this work we have presented an extended set of viscous hydrodynamic equations, in 
which particle number changing processes are taken into account via a rate equation. The 
rate equation of the particle number density is derived from the Boltzmann equation employ- 
ing the Grad's approximation. A one-dimensional boost-invariant expansion is considered 



17 





1.6 


"D 




CO 


1 h 


O 




P 




"D 


1 4 


h- 




CL 




■a 




h- 


1.3 


Q. 




7" 




T3 


1.2 


&3 




D_ 




^ 


1 1 


< 




CO 




!=■ 


1 


X3 


H 




Q. 




■a 


09 


h- 




Q. 




-z. 


0.8 



0.7 



x = 0.41 fm/c « x n - --- — 

x=0.8 fm/c = 2 x 

x=2.0 fm/c = 5 x - 




VI 




Jtf 


"1 


II 




[fllM 



0.5 



1.5 



2 2.5 

Pt (GeV) 



3.5 



4.5 



FIG. 11: (Color online) Ratio of transverse momentum spectrum extracted in BAMPS to the one 
calculated using the Grad's approximation at different times with r\j s = 0.35. 

for simplifying the derivation of the hydrodynamic equations. We have demonstrated that 
a proper treatment of particle number evolution is essential for a hydrodynamic description 
of the gluon system, especially for describing chemical equilibration and for determining the 
temperature and the momentum spectra of the system. 

We have compared the results between the hydrodynamic and transport calculations 
including inelastic 2 -H- 3 processes. The energy density and the shear pressure obtained 
from the third-order hydrodynamic equations agree well with the results from the transport 
approach using BAMPS even at large rj/s = 0.75, while the results from the Israel- Stewart 
hydrodynamics deviate from the BAMPS results by 10% — 20% for rj/s = 0.35 — 0.75. Both 
the IS and the third-order hydrodynamic calculations fail to meet the BAMPS results on the 
particle number density at large rj/s. The reason is that at large rj/s the viscous effect is so 
large that the Grad's ansatz for the off-equilibrium distribution should be modified beyond 
the second order in momentum, in order to give a more accurate rate equation of the particle 
number density. This is a highly interesting and timely subject for future investigations. 

The parton cascade BAMPS is presently being extended by including light quark degrees 
of freedom 33]. Light quarks are expected to equilibrate both kinetically and chemically 
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slower than a pure gluon system. Kinetic and chemical off-equilibrium effects will have 
significant influence on the jet and bulk medium physics. They will also be highly relevant 
concerning the slopes of transverse spectra and the yields of dileptons and photons, which 



were measured at RHIC 



34| to extract the initial temperature for hydrodynamic approaches. 
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